Universality of the edge tunneling exponent of fractional quantum Hall liquids 
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Recent calculations of the edge tunneling exponents in quantum Hall states appear to contradict 
their topological nature. We revisit this issue and find no fundamental discrepancies. In a micro- 
scopic model of fractional quantum Hall liquids with electron-electron interaction and confinement, 
we calculate the edge Green's function via exact diagonalization. Our results for v = 1/3 and 2/3 
suggest that in the presence of Coulomb interaction, the sharpness of the edge and the strength 
of the edge confining potential, which can lead to edge reconstruction, are the parameters that are 
relevant to the universality of the electron tunneling I-V exponent. 



One of the most intriguing characteristics of incom- 
pressible quantum Hall fluids is the nature of their edge 
excitations. Wenjj has argued that Hall fluids, which 
have no order parameter associated with a broken sym- 
metry, are nonetheless ordered topologically. While di- 
rect experimental probing of the topological order is dif- 
ficult, an indirect probe is provided by tunneling into the 
edge of the Hall fluid. In fact, by virtue of the topolog- 
ical order, edge modes are uniquely determined by the 
physics of the bulk and, in Abelian Hall states, form 
chiral Luttinger liquids 2] (CLL). For tunneling from a 
3-d Fermi liquid into the edge, CLL theory leads to a 
non-Ohmic tunneling current voltage relation / ~ V" , in 
sharp contrast to the Ohmic prediction of a Fermi-liquid 
dominated edge. 

For the Hall states dX v = n/{2np + 1) (where n is a 
nonzero integer and p is an even positive integer), the 
edge for n > does not contain counter propagating 
modes and the exponent is a = p + 1, independent of 
n. The situation is more complicated for n < where 
counter propagating modes can be back-scattered. How- 
ever, in the presence of disorder, the exponent takes on 
the universajjal value a = p + 1 — 2/|n|. While experi- 
ments 0, 0, Us ISl have confirmed the nontrivial power- 
law behavior, they do not agree with CLL values 2]. In 
particular, one experiment found an approximate de- 
pendence oi aK, 1/u. 

Earlier attempts U H [H E [13 to resolve the ap- 
parent discrepancy between experiment and theory have 
been summarized in Many of these approaches 

have invoked additional physics beyond the standard the- 
ory to address the shortcoming rather than invalidating 
the basic CLL picture. One such addition arises from 
the presence of a positive background charge. On purely 
electrostatic grounds the electron density near the edge 
may become quite different from that of an ideal edge [a| ■ 
This effect can even lead to the reconstruction of the 
edge 0, 0, 0, [3 provided the background charge is 
sufficiently far from the electron layer (which is usually 
the case in cleaved-edge samples). As a consequence, the 
tunneling characteristics could become very sensitive to 
the edge profile and the universal tunneling characteris- 



tics may not necessarily be observed. 

Meanwhile, yet another line of thought 0, 113, 0| 
which questions the role of the range of electron-electron 
interaction has emerged. Tsiper and Goldman (TG) 
studied the edge wave function using exact diagonaliza- 
tion in the presence of Coulomb interaction [ijj. They 
concluded that the tunneling exponent depends on the 
range of interaction. Crucial to their study is the as- 
sumption that the exponent a may be obtained from the 
ratio of the electron occupation numbers of the two out- 
ermost occupied orbitals for the corresponding Laughlin 
state in the disk geometry, i.e. a = p{m^^^—l)/p{m^^^^). 
This relation, however, has been derived only in the case 
of ultra short range interactions and its validity for the 
more generic finite range case has not been established. 

U sing composite fermion (CF) theoryji^, Mandal and 
Jain[2^ |2J (MJ) have arrived at essentially the same 
conclusion. These authors adopted a hard edge by cut- 
ting off angular momentum larger than rrimax = 3(Af — 1) 
for = 1/3 and, as TG, ignored the background charge. 
They found that for the ultra short range potential 
(which produces the Laughlin state), the asymptotic edge 
Green's function exponent agreed with CLL theory. On 
the other hand, for generic potentials, in particular the 
Coulomb potential, a substantial reduction of the expo- 
nent from the CLL value of 3 was observed. MJ at- 
tributed this reduction to the residual repulsion (beyond 
their hard core) among the composite fermions generated 
by the long-range Coulomb potential. For v = 1/3, the 
exponent is below 2.5 and even larger reductions were 
found for ly = 2/5 and 3/7. 

These results not only are at odds with the predictions 
of CLL, but cast doubt on the most crucial element of 
the FQH physics itself, namely the concept of topological 
order. The unusual properties of the chiral edge liquid is 
understood to be the signature of the topological struc- 
ture of the bulk and therefore should persist as long as 
the bulk exhibits the FQH effect 0. Hence, one expects 
the same exponent irrespective of the range of the inter- 
actions so long as the bulk physics remains the same. 

In this paper we show that there are no fundamen- 
tal contradictions with CLL and/or the topological order 
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of FQH states. In the presence of long-range Coulomb 
interaction, our findings suggest that the detail of the 
edge confinement is relevant to understanding the be- 
havior of the edge tunneling exponent. We first address 
the edge exponent in a system with long-range Coulomb 
interactions in the absence of neutralizing background 
charge. To this end we evaluate the edge Green's func- 
tion by exact diagonalization in a microscopic model of 
the FQH liquids. We impose an edge confining potential 
by restricting the single-particle angular momenta to be 
< JTimax- We find that, for v = 1/3, the tunneHng ex- 
ponent remains unchanged with Coulomb interaction for 
soft edges (large mmax)- This is in sharp contrast to the 
reduction of a as found previously by MJ for hard edge 
confinement (small mmax)- We then investigate the effect 
of the edge potential induced by background charge in the 
presence of long-range interaction. For = 1/3 and weak 
confining potential, we again observe substantial devia- 
tions from the universal value for hard edges, which rnay 
be highly relevant to the experimental studies |Sl la LD| ■ 
We also find finite-size corrections to a for soft edges, con- 
sistent with the edge reconstruction scenario jlSl Il6l llq . 
For V = 2/3, we find behavior consistent with strongly 
coupled edges for strong confining potential and with an 
uncoupled outer v = \ edge for weak confining poten- 
tial. Finally, we further emphasize the importance of 
including the edge confining potential by comparing the 
ground state quantum numbers for z/ — 2/5 obtained by 
exact diagonalization to those of the CF and hierarchical 
constructions. 

We consider a microscopic model of a two-dimensional 
electron gas (2DEG) confined to a two-dimensional disk 
with neutralizing background charge distributed uni- 
formly on a parallel disk of radius a, at a distance d 
above the 2DEG. The choice of a = ^j2N/v guaran- 
tees that the disk encloses N electrons and exactly N/v 
magnetic flux quanta for the desired filling factor v. 
The bare Coulomb interaction between the background 
charge and the electrons gives rise to the confining po- 
tential. We use the same Hamiltonian as in our previous 
studvp^. We confine the electrons to the lowest Landau 
level (LL) and employ symmetric gauge wavefunctions: 
(/)„(z) = (27r2™TO!)-i/2z™e-l^l'/4^ where z = x + iy \s 
the complex coordinate. In this paper the distances are 
measured in units of the magnetic length Ib = ^Jh/eB. 

We diagonalize the Hamiltonian to obtain the exact 
many-body ground state if) using the Lanczos algorithm. 
We then calculate the equal-time edge Green's function, 



(1) 



where vE'Ki") ^^'^ ^e(r') are field operators, which create 
and annihilate an electron at r and r', respectively, on the 
edge of the 2DEG disk with a radius of R and |r — r'| = 
2Rsm{9/2). The choice of R is not crucial and will be 




0.6 0.7 0.8 0.9 1 
|sin(e/2)| 





0.5 0.6 0.7 0.8 0.9 1 
|sin(e/2)| 



FIG. 1: The edge Green's function \G\ for the Laughlin state 
with 6-9 electrons at filling fraction u — 1/3 (a) as a function 
of 6 and (b) as a function of |r — r'|. (c) |G'(S)| for 8 electrons 
with Coulomb interaction confined to orbitals with the largest 
angular momentum mmax = 23-29. (d) |G(|r — r'j)| for A*" = 
6-9 electrons with Coulomb interaction confined to orbitals 
with rriniax ~ {N — + 5 for u = 1/3. The lines in the 
log- log plots (a)-(d) correspond to a power-law behavior with 
a = 3.2. 



specified later. In the large |r — r'| limit, the edge Green's 
function is expected to exhibit the asymptotic behavior 



\Gedge{r - r' 



cx |sin(6'/2)|- 



(2) 



Because of the relativistic invariance of CLL, the equal 
time and equal distance exponents of the Green's func- 
tion are equal; the latter is measured in tunneling exper- 
iments. 

For comparison, we first consider the ultra short-range 
hardcore potential, for which the Laughlin state is the 
exact ground state. We do not include the background 
confining potential, but choose the ground state with the 
appropriate total angular momentum. Figure ^a) shows 
the edge Green's function (i? = y^2N/i/) for the Laughlin 
state with 6-9 electrons at filling fraction i' — 1/3. We 
use least-square fit to match our data to the power-law 
\G{e)\ - |sin(6l/2)|-" close to lsin(6l/2)l = 1, and obtain 
a — 3.2 ± 0.2. The errorbar reflects the dependence of 
a on system size and range of data to fit. This result 
is in good agreement with a = 3 as predicted by the 
CLL theory. |G'(6')| for = 6 shows weak oscillation 
around the power-law fitting curve, but the finite-size 
effects become very weak for N > 7. In Fig. Cfb), we 
replot \G\ as a function of |r — r'|. We observe perfect 
scaling even for distances |r — r'| as small as one magnetic 
length, which is a strong indication that finite size effects 
are indeed negligible. 

Next we consider the long-range Coulomb interaction. 
There is an important difference here with the Laughlin 
state so far as the edge is concerned; in the latter there 
are no occupied single-particle angular momenta that ex- 
ceed TO^ax = (-^ ~ l)/'^- Thus we need to enlarge our 
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basis set and find the number of orbitals beyond which 
the properties of the system converge. Figure Q^c) plots 
\Gedge\ for the Coulomb interaction and N — 8 electrons 
at filling fraction = 1/3 for an increasing number of 
orbitals (mmax + 1 since we label from m = 0). We de- 
fine the edge by choosing R ~ ^2(m„iax + 1) hereafter. 
For TTimax < 26 (hard edge), we find a weak oscillation 
of |Gedge| cvcn near the largest distance of the system. 
These oscillations are probably induced by the compe- 
tition between the long-range interaction and the edge 
confining potential. Similar oscillations, existing gener- 
ically at other filling fractions, can also be observed in 
the electron density profile in the presence of Coulomb 
interaction Therefore, fitting \Gedge\ to Eq. (O 

to extract a may not produce an accurate exponent. On 
the other hand, for mmax > 26 (soft edge), jGedgel can be 
fit very well by the power-law with a — 3.2 ±0.2, which is 
the same as the ultra short-range interaction exponent. 
In Fig. njd), we again show a scaling plot of \G\ over 
|r — r'l for = 6 — 9 electrons and Wmax = {N — l)/i' + 5 
at = 1/3, again for the Coulomb interaction. Even with 
long range interaction, the data shows good scaling with 
only small deviations at length scales below 8. We note 
that the choice of Wmax here is the same as in Ref. [iflj . 
The difference in the exponent is caused by the manner it 
was determined. We have verified that the formula used 
by TG does not agree with the exponent in the Green's 
function. 

So far we have excluded the background confining po- 
tential. Without the background charge, electrons tend 
to move to the edge to reduce their Coulomb repulsion. 
This seems to induce strong density oscillations near 
the edge, extending into the bulk rather than forming a 
roughly uniform droplet. Nor does it conform to the ex- 
periments where a confining potential is always present. 
Figure |2a) shows the edge Green's function for 8 elec- 
trons with mmax = 23 (hard edge) with the correspond- 
ing confining potential for v = 1/3. For d = 1.0, where 
there is no edge reconstruction (strong confining poten- 
tial), we find that G{9) agrees very well with a power 
law of G ~ |sin(6'/2)|" with a = 3.2 ± 0.1. This is 
equal to the exponent in the complete absence of any 
confining potential. This is because the background has 
largely mitigated the combined effects of the long-range 
repulsion and the hard-edge confinement. However, for 
d > dc « 1.5, due to edge reconstruction, this is no 
longer the case. For d = 1.8 > dc (weak confining po- 
tential), G{d) increases its value as a result of electrons 
moving closer to the edge and changes a to 2.2 ± 0.1. 
Again, one can see this qualitatively on the electrostatic 
level; the electron occupation ratio in the lowest Lan- 
dau level near the edge is larger than 1/3. We next re- 
lax the cutoff in angular momentum space and compare 
G{e) with two different m^ax for d = 1.8 (Fig. E^b)). 
For mmax = 29 (soft edge), g{6) shows a crossover from 
a power law with a ~ 5.0 to one with a ~ 3.0 near 
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FIG. 2: The edge Green's function \G\ of 8 electrons at 
1/ = 1/3 with Coulomb interaction and the background charge 
confining potential for (a) mmax ~ 23 before (d — 1.0) and 
after the edge reconstruction (d = 1.8), and (b) rrimax ~ 23 
and 29 with d = 1.8. The straight lines are power- law fits 
with a — 2.2 and 3.2 in (a), and a — 2.2 and 3.0 in (b), 
respectively. 

sin(0/2)| « 0.75. This suggests that the true asymptotic 
behavior in the reconstructed case can only be observed 
at a larger length scale. Such behavior agrees qualita- 
tively with the edge reconstruction corrections to a at 
short distances: Sa oc v^/v^ [IE 113 j where and v are 
velocities of neutral and charge modes, respectively. 

The significant drop in a for v = 1/3 in the case of hard 
edges corroborates the previous results for long-range in- 
teraction Efi E]| . While it is impossible to deter- 
mine with certainty what happens in the thermodynamic 
limit, we agree with the assessment of MJ that these re- 
ductions are not finite-size artifacts, notwithstanding the 
large distance oscillations we find in G. However, our re- 
sults for soft edges appear to show that the non-universal 
behavior has more to do with the details of edge confine- 
ment than the range of the interaction potential. Indeed, 
Figs. Hand 13 suggest that the edge confinement, through 
"T-max as well as d, is relevant to a in the presence of 
long-range interaction. These issues are moot for the ul- 
tra short-range interaction (unless mmax < {N ~ l)/i^, in 
which case the Laughlin state cannot even be realized). 
As pointed out in Ref. the CF ground state with 
one CF exciton involves only single particle states with 
m < mmax = 3(7V — 1), corresponding to the hard edge 
in our study. It would be interesting to find the precise 
CF state that would correspond to our soft edge profile. 

We have also studied the behavior of the edge Green's 
function at other filling fractions, such as u — 2/3, 
which is not investigated in the CF approach of MJ. The 

= 2/3 droplet can be regarded as a = 1/3 hole droplet 
superimposed on a u = 1 electron droplet. It therefore 
supports an inner v — t/i edge and an outer v ^ 1 
edge Figure 131a) compares the edge Green's func- 

tion for 18 electrons in 27 orbitals (hard edge) with the 
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FIG. 3: (a) The edge Green's function \G\ for u = 2/3 with 
N = 18 and rrimax = 27. The straight lines are power-law fits 
with exponent a = 1.4 and 1.0 for d — 0.2 and 2.0, respec- 
tively, (b) The total angular momentum Mo of the ground 
state as a function of d, the distance between charge layers 
for A*' = 10 electrons at v — 2/5 (mmax = 25). 



Ref. 28]). Evidently, the CPs in our case either form 
partially filled Landau levels and/or have reconstructed 
edges. We point out that the ground state with Mq = 135 
near d — 2.5 is a = 1/3 Laughlin state (a filled lowest 
Landau level of CFs). The transition from i' — 2/5 to 
v — 1/3 as d increases is an artifact of our finite system 
size. Therefore, the further increase of Mq for even larger 
d is, in fact, the consequence of the reconstruction of the 
I' = 1/3 edge 0,^^, albeit we set up the background 
charge distribution for v = 2/5. We note that the numer- 
ical approach applies to arbitrary u, which may involve 
many CF landau levels or be nested deep in the hierarchy. 

We thank Matt Grayson, Jainendra Jain, Akakii Me- 
likidze, Xiao-Gang Wen and Kun Yang for helpful dis- 
cussions. The work is supported by the Schwerpunk- 
tprogramme "Quanten-Hall-Systeme" der DFG and by 
US DOE under contract DE-FG03-02ER-45981. 



corresponding confining potential for d = 0.2 and 2.0. 
For strong confining potential {d — 0.2), we find, by fit- 
ting G{d) to a power law, that a — lA regardless of 
TOmax- This is close to 1/v — 1.5 and we speculate that 
the two counter propagating edge modes strongly couple 
and reconstruct into a dominant charge mode and a neg- 
ligible neutral mode. On the other hand, for weak confin- 
ing potential {d = 2.0), we find a = 1.0, which probably 
is the fingerprint of the reconstructed outer edge of the 
y = 1 fluid. 

We have demonstrated the nontrivial effects of the edge 
confining potential. We would like to emphasize that 
the inclusion of the realistic confining potential not only 
guarantees the charge neutrality and the homogeneity of 
the 2DEG, but also provides a numerical method of de- 
termining the total angular momentum of the most sta- 
ble states. The ground-state angular momentum Mq of 
the interacting system is known for the principal filling 
fractions, such as v = 1/3, for which Laughlin's varia- 
tional wave function |25j is a good approximation. Based 
on Haldane's hierarchical construction [2^ or Jain's CF 
theory (2 the variational wave function for certain 
filling fractions such as = 2/5 can also be written 
down [Ullil. In Fig. ^h), we plot Mq as a function 
of d, the distance between charge layers, for a v = 2/5 
quantum Hall liquid with TV = 10 electrons. Here, Mq is 
determined by looking for the lowest ground state energy 
in all angular momentum subspaces, since the system 
maintains rotational symmetry. The nondecreasiiig curve 
is similar to those found at other filling fractions 
We note that Mq takes a discrete value of 111, 119, 127, 
and 135 for d = 0.2-3.0, while the CF or the hierarchy 
construction predicts Mq = 111, 117, 125, and 135 for 3, 
2, 1, and CFs in the second CF Landau level (forming 
maximum density droplets), respectively |2^ (see also 
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